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I. INTRODUCTION 



An analytical solution of the radial Schrodinger equation is of high importance in nonrel- 
ativistic quantum mechanics, because the wave function contains all necessary information 
for full description of a quantum system. There are only few potentials for which the radial 
Schrodinger equation can be solved explicitly for all n and /. So far, many methods were 
developed, such as supersymmetry (SUSY) [1,2] and the Pekeris approximation [3-6], to 
solve the radial Schrodinger equation exactly or quasi-exactly for Z 7^ within these poten- 
tials. Levai and Williams suggested a simple method for constructing potentials for which 
the Schrodinger equation can be solved exactly in terms of special functions [7] and showed 
relationship between the introduced formalism and supersymmetric quantum mechanics [1]. 
Amore et al. presented a new method [8] for the solution of the Schrodinger equation appli- 
cable to problems of nonperturbative nature. In addition, they applied the method to the 
quantum anharmonic oscillator and found energy eigenvalues and wave functions, even for 
strong couplings. 

The radial Schrodinger equation for the Woods-Saxon potential [9] can not be solved 
exactly for / 7^ 0. However, Fliigge gave an exact expression for the wave function, but a 
graphical method was suggested for the energy eigenvalues at / = [4]. It is well known 
that the Woods-Saxon potential is one of the important short-range potentials in physics. 
Furthermore, this potential were applied to numerous problems, in nuclear and particle 
physics, atomic physics, condensed matter, and chemical physics. 

Recently, an alternative method known as the Nikiforov-Uvarov (NU) method [10] was 
proposed for solving the Schrodinger equation. Therefore, it would be interesting and im- 
portant to solve the nonrelativistic radial Schrodinger equation for Woods-Saxon potential 
for I 7^ 0, since it has been extensively used to describe the bound and continuum states of 
the interacting systems. Thus, one needs to obtain the energy eigenvalues and correspond- 
ing eigenfunctions of the one particle problem within this potential. The NU method was 
used by C. Berkdemir et al. [11] to solve the radial Schrodinger equation for the general- 
ized Woods-Saxon potential for I = 0. However, it this work, the authors made errors in 
application of the NU method, which led to no-correct results [12]. 

In this work, we solve the radial Schrodinger equation for the standard Woods-Saxon 
potential using NU method [10], and obtain the energy eigenvalues and corresponding eigen- 
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functions for arbitrary / states. 



II. NIKIFOROV-UVAROV METHOD 

The Nikiforov-Uvarov (NU) method is based on the solutions of general second-order 
linear equations with special orthogonal functions. It has been extensively used to solve 
the nonrelativistic Schrodinger equation and other Schrodinger-like equations. The one- 
dimensional Schrodinger equation or similar second-order differential equations can be writ- 
ten using NU method in the following form: 

Wz) + ^(*) + Hl^) = 0, (2.1) 

where a(z) and a(z) are polynomials, at most second-degree, and t(z) is a first-degree 
polynomial. 

Using Eq.(2.1) the transformation 

,/>(z) = Hz)y(z) (2.2) 
one reduces it to the hypergeometric-type equation 

a(z)y" + r(z)y' + Xy = 0. (2.3) 



The function <&(z) is defined as the logarithmic derivative [10] 

&(z) 7T(z) 



(2.4) 



$(z) <j{zY 

where ■k(z) is at most the first-degree polynomial. 

The another part of ip(z), namely y(z), is the hypergeometric-type function, that for fixed 
n is given by the Rodriguez relation: 

where B n is the normalization constant and the weight function p(z) must satisfy the con- 
dition [10] 

±(a(z)p(z))=T(z)p(z), (2.6) 

with t(z) = t(z) + 2n(z). 
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For accomplishment of the conditions imposed on function p(z), the classical orthogonal 
polynomials, it is necessary, that polynomial t(z) becomes equal to zero in some point of an 
interval (a, b) and derivative of this polynomial for this interval at a(z) > will be negative, 
i.e. r'(z) < 0. 

The function tt(z) and the parameter A required for this method are defined as follows: 



n(z) = ± J (^y^) '-v + k*, (2.7) 

\ = k + n\z). (2.8) 

On the other hand, in order to find the value of k, the expression under the square root 

must be the square of a polynomial. This is possible only if its discriminant is zero. Thus, 
the new eigenvalue equation for the Schrodinger equation becomes [10] 

A = A n = -nr' - n{n - l \ »^ ( n = 0, 1, 2, ...). (2.9) 
After the comparison of Eq.(2.8) with Eq.(2.9), we obtain the energy eigenvalues. 



III. SOLUTIONS OF THE SCHRODINGER EQUATION WITH THE WOODS- 
SAXON POTENTIAL 

The standard Woods-Saxon potential [9] is defined by 

V(r) = - V ° (r . RoV a « Ro- (3.1) 
1 + exp y 1 —^) 

This potential was used for description of interaction of a neutron with a heavy nucleus. The 
parameter Rq is interpreted as radius of a nucleus, the parameter a characterizes thickness 
of the superficial layer inside, which the potential falls from value V = outside of a nucleus 
up to value V = —Vq inside a nucleus. At a = 0, one gets the simple potential well with 
jump of potential on the surface of a nucleus. 

The radial Schrodinger equation [13, 14] with Woods-Saxon potential is 



d 2 R 2 dR 2a 
\ h — 

dr 2 r dr h 2 



1 + exp (r=f*) 



R- ^j-^ fl = 0, (0<r <oo), (3.2) 



where I is the angular momentum quantum number and fi is the reduced mass. 
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Introducing a new function 

u(r) = rR(r), 

Eq.(3.2) takes the form 



d 2 u(r) 2/i 
dr 2 ft 2 



v h 2 i{i + \) 

hi + 



u(r) = 0. (3.3) 



1 + exp ( r ^ a ) 2//r 2 
Equation (3.3) has the same form as the equation for a particle in one dimension, except for 
two important differences. First, there is a repulsive effective potential proportional to the 
eigenvalue of h 2 l(l + 1). Second, the radial function must satisfy the boundary conditions 
«(0) = and u(oo) = 0. 

It is sometimes convenient to define in Eq.(3.3) the effective potential in the form: 



VefAr) = V(r) + ^±^-. (3.4) 



Then, the radial Schrodinger equation given by Eq.(3.3) takes the form 

^ + %[E-V eff (r)]u = 0. (3.5) 

If in Eq.(3.1) introduce the notations 

_ r - Rq _ Ro 

X , OL , 

Rq a 

then the Woods-Saxon potential is given by the expression 

V - Vo 

vws - --r- 7 — r • 

1 + exp (ax) 

The effective potential together with the WS potential for I ^ can be written as 



V ef f(r) = V t (r) + V ws (r) = ^f±R _ V ° (3.6) 

2^r 2 1 + exp(ax) 



It is known that the Schrodinger equation cannot be solved exactly for this potential at the 
value / 7^ using the standard methods as SUSY and NU. From Eq.(3.6) it is seen that the 
effective potential is a combination of the exponential and inverse square potentials, which 
cannot be solved analytically. Therefore, in order to solve this problem we can take the most 
widely used and convenient for our purposes Pekeris approximation. This approximation is 
based on the expansion of the centrifugal barrier in a series of exponentials depending on the 
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internuclear distance, taking into account terms up to second order, so that the effective l- 
dependent potential preserves the original form. It should be pointed out, however, that this 
approximation is valid only for low vibrational energy states. By changing the coordinates 



x 



r-Rp 



or r = Rq(1 +x), the centrifugal potential is expanded in the Taylor series around 



the point x = (r = Rq) 

h 2 l(l + 1) h 2 l(l + 1) 1 



V l (r) = 



where 6 = 



2^r 2 



2^R 2 (1 + xf 



5(l-2x + 3x 2 -4x 3 + ...) , (3.7) 



According to the Pekeris approximation, we shall replace potential V/(r) with expression 



vr(o = *(co+ r 



Ci 



Co 



(3.8) 

+ expax (1 + expax) 
In order to define the constants Co, C\ and C 2 , we also expand this potential in the Taylor 
series around the point x = (r = Rq): 



„ Ci C 2 \ Q! . a 2 n a 3 , a 4 4 

Co + y + T J - j (Ci + C 2 ) x + — C 2 x 2 + — (d + C 2 ) x 3 - — C 2 x 4 + 



(3.9) 

Comparing equal powers of x Eqs.(3.7) and (3.9), we obtain the constants C , C\ and C 2 : 

4 12 8 48 48 

<-^0 — 1 1 nl W — ^) ^2 — ~ ^- 

a cr a ar cr 

Now, the effective potential after Pekeris approximation becomes equal to 

v e * ff (x) = v;{x) + v W3 {x) = sc - : v ° - 6Cl ; + 



(3.10) 



l + exp(ax) ' (1 + exp(ax)) 2 ' 

Instead of solving the radial Schrodinger equation for the effective Woods-Saxon potential 
V e ff(r) given by Eq.(3.6), we now solve the radial Schrodinger equation for the new effec- 
tive potential K//( r ) gi ven by Eq.(3.9) obtained using the Pekeris approximation. Having 
inserted this new effective potential into Eq.(3.5), we obtain 



d 2 u 2/i 
dr 2 h 2 



E-SC + 



Vo - 8d 



5C, 



l + exp(^) (i + ex p(^)) : 
We use the following dimensionless notations: 



u = 0. 



(3.11) 



2/i (E - SCq) a 2 

h 2 



; P 2 = 



2fi (Vo - 6d) a 2 
h 2 



7 = 



2^C 2 a 2 



(3.12) 
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with real e > for bound states; (3 and 7 are real and positive. 
If we rewrite Eq.(3.11) by using a new variable of the form 

r-RoW 1 



z = yl + exp 
we obtain 

1 - 2z -e 2 + B 2 z - Vz 2 

u'\z) + -t^—M + +P J u(z) = 0,(0<z< 1), (3.13) 
z{l-z) [z(l-z)) z 

with t(z) = 1-2z; a(z) = z(l - z); a(z) = -e 2 + (3 2 z - 7V. 
In the NU-method, the new function ir(z) is 

n(z) = ±v/e 2 + (k - (3 2 ) z - {k - 7 2 ) z 2 . (3.14) 

The constant parameter k can be found employing the condition that the expression 
under the square root has a double zero, i.e., its discriminant is equal to zero. Hence, there 
are two possible functions for each k: 

( (e - Je 2 -(3 2 + 7 2 ) z-e, fork = (3 2 - 2e 2 + 2eJt 2 - f3 2 + 7 2 , 

n(z) = ±{ } V , J V , (3.15) 

'[ U + Ve 2 - (3 2 + 7 2 J z-e, for k = [3 2 - 2e 2 - 2e^e 2 - (3 2 + 7 2 . 

According to the NU-method, from the four possible forms of the polynomial 7r(z) we select 
the one for which the function r(z) has the negative derivative and root lies in the interval 
(0,1). Therefore, the appropriate functions ir(z) and t(z) have the following forms: 

n(z) = e - (e + ^e 2 - /3 2 + 7 2 ) z, (3.16) 

t(z) = 1 + 2e - 2 + e + ^e 2 - f3 2 + 7 2 ) z, (3.17) 

and 

k = f3 2 - 2e 2 - 2e v / e 2 - (3 2 + 7 2 . (3.18) 
Then, the constant A = A; + 7t'(2;) is written as 



A = (3 2 - 2e 2 - 2e v / e 2 - f3 2 + ~/ 2 - e + ^e 2 - f3 2 + 7 2 . (3.19) 
An alternative definition of A n (Eq.(2.9)) is 

A = A n = 2 (e + x/e 2 - (3 2 + 7 2 ) n + n{n + 1). (3.20) 
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Having compared Eqs.(3.19) and (3.20) 

/5 2 -2e 2 -2e v / e 2 - (3 2 + 7 2 -e- ^e 2 - (3 2 + 7 2 = 2 (e + v^ 2 ~ /? 2 + 7 2 ) ra+ra(ra+l), (3.21) 
we obtain 

e + A / P 2 _ ^2 4^2 4. „_ + 

2 2 



e + Ve 2 - P 2 + 7 2 + n + \ - \ v 7 ! + 4 7 2 = 0, (3.22) 



or 



e + Ve 2 - /3 2 + 7 2 - n' = 0. (3.23) 

Here 

n ,_ n+ yr±^l (3 . 24) 

n being the radial quantum number (n = 0, 1, 2, . . .). From Eq.(3.23), we find 
Because for the bound states e > 0, we get 

ri > (3.26) 

and 

-n' 2 </? 2 - 7 2 <n' 2 . (3.27) 

If n' > 0, there exist bound states, otherwise, there are no bound states at all. By using 
Eq.(3.24) this relation can be recast into the form 



Jl + 4 7 2 - 1 

< n < -¥ , (3.28) 

i.e. it gives the finite coupling value. 

If — n' 2 < [3 2 — 7 2 < n' 2 , there exists bound states; otherwise there are no bound states. 
Inequality, which obtained after substituting (3, 7, Co, Ci, C2 into Eq.(3.27), gives the definite 
coupling value for the potential depth Vq\ 

M Ji—n' 2 < V < — -^—n' 2 (3 29) 

a 2fj,a 2 a 2/ia 2 

After substituting a,7,5,C 2 and Eq.(3.24) into Eqs.(3.28),(3.29), we find 



^ . 192a 4 /(Z+l) -1 
V i + K? 1 



< n < (3.30) 



and 



4h 2 al{l + 1) _ ( \ | 192a 4 /(/ + l) 



/iRq 8^a 2 



i?0 



2n - 1 < 



or 



T , 4n 2 a/(Z + l) ft 2 / / 192a 4 /(/ + l) n 



4^ 2 Z(Z + 1) a £ 2 R 2 ( I 192aH(l + 1) 



2n - 1 < 



Vo^< 



4£ 2 Z(Z + 1) a n 2 i? 2 / f 192a 4 /(/ + 1) 

— s — i^ + 8^v — ^ — 2n -^ 



By defining 



(3.31) 



(3.32) 



Omin 



Ah 2 al(l + 1) /l 2 



hRq 8/ia 2 



, i + CT + i)_ 2n _ i 

/t 



(3.33) 



and 



Vornax 



4fr 2 a/(Z + 1) ft 2 



' 8/ia 2 

Eq.(3.31) takes the following form: 



192a 4 /(/ + 1) n 

L + 3 -"2n- I 

K 



(3.34) 



Vomin < Vq < Vo max 



(3.35) 



The exact energy eigenvalues of the Schrodinger equation with the Woods-Saxon potential 
are derived as 

/ n '2 + g2_ 2\2 

^ = fl?o - (V - 5C\) J , ) ■ (3.36) 



2(3n' J 

Substituting the values of 5, C , Ci, C 2 , n', /3 and 7 into (3.36), one can find E ni 

h 2 l(l + l) ( 12a 2 



2fiR 2 



Rl 



h 2 



1 + ^ - 2n - 1 



2/za 2 



16 



+ 



/ia 2 V 4Z(Z+l)a 3 
"IP 



Z? 3 

-rt 



1 + ""(ft 1 * 4 - 2n - 1 



+ 



[iVoa 2 



(3.37) 



If both conditions (3.28) and (3.29) ((3.30) and (3.31)) are satisfied simultaneously, the 
bound states exist. Thus, the energy spectrum equation (3.37) is limited, i.e. we have only 
the finite number of energy eigenvalues. 
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For very large V the /-dependent effective potential has the same form as the potential 
with / = 0. From Eqs.(3.30) and (3.31), it is seen that if I = 0, then one gets n < 0. 
Hence, the Schrodinger equation for the standard Woods-Saxon potential with zero angular 
momentum has no bound state. For larger values of Vq (Vq > Vo max ), the conditions (3.30) 
and (3.31) are not satisfied. Therefore, no bound states exist for larger values of Vq (Vq > 
Vomax)- Similarly, no bound states exist for smaller values of Vq (0 < Vq < Vomm) too. 

According to Eq.(3.37), the energy eigenvalues depend on the depth of the potential V , 
the width of the potential Rq, and surface thickness a. Any energy eigenvalue must be 
less than Vq. If constraints imposed on n and Vq satisfied, the bound states appear. From 
Eq.(3.31), it is seen that the potential depth increases when the parameter R increases, but 
the parameter a is declined and vice versa. Therefore, one can say that the bound states 
exist within this potential. 

By using the empirical values r = 1.285/m and a = 0.65/m taken from Ref.15, the 
potential depth V = (40.5 + 0.13A)MeV = 47.78MeV and the radius of the nucleus R = 
tqA 1 ^ = 4.9162/m are calculated for the atomic mass number of target nucleus A = 56 
and the reduced mass n = 0.50433-u. In Table 1, energies of the bound states obtained 
numerically for the spherical standard Woods-Saxon potential for some values of / and n 
are given. It is clearly seen from Table 1 that the Vo m i„ and E n \ increase, but VQ ma x reduces 
when for the fixed values / the parameter n is increased. Therefore, AV = V 0max — V 0m i n 
reduces, when n increases. 

In addition, we have seen that there are some restrictions on the potential parameters 
for the bound state solutions within the framework of quantum mechanics. That is, when 
the values of the parameters Vq and n satisfy the conditions (3.30) and (3.31), we ob- 
tain the bound states. We also point out that the exact results obtained for the standard 
Woods-Saxon potential may have some interesting applications for studying different quan- 
tum mechanical and nuclear scattering problems. Consequently, the found wave functions 
are physical ones. 

Now, we are going to determine the radial eigenfunctions of this potential. Having sub- 
stituted ty(z) and a(z) into Eq.(2.4) and then solving first-order differential equation, one 
can find the finite function $(z) in the interval [0, 1] 

= z € (l-z)^ e2 -^ 2 . (3.38) 
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It is easy to find the second part of the wave function from the definition of weight 
function: 

p(z) = z^{\-zf^^^\ (3.39) 



and substituting into Rodrigues relation (2.4), we get 

(3.40) 



y n {z) 



B n z~* (1 - z)- 2 ^™* 



dz n 



y n+2e ^ _ ^n+2 v 4 2 -/3 2 + 7 2 



where B n is the normalization constant and its value is ^ [16]. Then, y n is given by the 
Jacobi polynomials: 



(2e,2 A /e 2 -/3 2 +7 2 ) 

y n (z) = p} v J (l-2z), 

where 

1 rl n v 

P<?'®(1 - 2z) = -z~ a (1 - z)"" -f- ^ n+a (1 - ^) n+/3 
The corresponding u n i(z) radial wave functions are found to be 

u nl {z) = C nl z^ (1 - z)^^P^^^\l - 2z), (3.41) 

where C n \ is a new normalization constant determined using f^°[u n i(r)] 2 dr = 1. 

IV. CONCLUSION 

In this paper, we have analytically calculated energy eigenvalues of the bound states 
and corresponding eigenfunctions in the new exactly solvable Woods-Saxon potential. The 
energy eigenvalue expression for Woods-Saxon potentials is given by Eq.(3.37). As it should 
be expected (see Eq.(3.37)), for any given set of parameters V ,Ro and a, the energy levels 
of standard Woods-Saxon potential are positive. We can conclude that our results are 
interesting not only for pure theoretical physicist but also for experimental physicist, because 
the results are exact and more general. 
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1 


n 


V 0min ,MeV 


Vomaxi MeV 


Vo, MeV 


E nh MeV 


1 





3.5590 


3.7191 


3.6 


2.3374 


2 





10.2654 


11.5690 


10.5 


7.0068 


3 





19.5398 


24.1289 


20 


14.0327 


4 





30.8571 


41.9241 


35 


22.6509 


5 





43.8248 


65.3471 


47.78 


34.7761 


6 





58.1722 


94.6684 


90 


35.3171 


7 





73.7179 


130.0696 


120 


46.5875 


8 





90.3395 


171.6730 


160 


54.5035 


9 





107.9536 


219.5620 


200 


67.4620 


10 





126.5020 


273.7949 


240 


85.1164 


11 





145.9433 


334.4129 


270 


102.6090 


12 





166.2474 


401.4463 


284 


153.2491 


12 


1 


282.9924 


284.7013 


284 


182.4285 


13 





187.3992 


474.9172 


330 


177.8142 


13 


1 


326.8691 


335.4402 


330 


212.6080 


14 





209.3607 


554.8423 


390 


198.7520 


14 


1 


371.7128 


392.4902 


390 


237.9290 


15 





232.1402 


641.2347 


450 


223.1066 


15 




417.4850 


455.8899 


450 


267.3521 


20 





357.9160 


1170.4901 


764.2 


390.3832 


20 


1 


659.4146 


868.9915 


764.2 


465.7578 


20 


2 


764.1028 


764.3034 


764.2 


491.9299 


30 





667.3430 


2716.9849 


1690 


834.2010 


30 


1 


1204.0683 


2180.2595 


1690 


968.3811 


30 


2 


1543.9833 


1840.3445 


1690 


1053.3543 


30 


3 


1687.0879 


1697.2400 


1690 


1088.9078 
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1 


n 


V 0min ,MeV 


V 0max ,MeV 


V ,MeV 


E nh MeV 


40 





1052.7566 


4915.3065 


3000 


1430.1249 


40 


1 


1826.2400 


4144.8220 


3000 


1623.4737 


40 


2 


2402.9130 


3565.1490 


3000 


1767.5873 


40 


3 


2782.7757 


3185.2864 


3000 


1862.3459 


40 


4 


2965.8279 


3002.2441 


3000 


1904.9234 


50 





1513.7067 


7765.9020 


4600 


2225.1111 


50 


1 


2524.5792 


6755.0294 


4600 


2477.7687 


50 


2 


3338.6413 


5940.9673 


4600 


2681.1670 


50 


3 


3955.8930 


5323.7156 


4600 


2835.2052 


50 


4 


4376.3344 


4903.2743 


4600 


2939.3913 


50 


5 


4599.9653 


4679.6434 


4600 


2986.8600 


100 





4948.2206 


31806.3078 


18400 


8461.6670 


100 


1 


7148.9350 


29605.5934 


18400 


9011.8467 


100 


2 


9152.8390 


27601.6894 


18400 


9512.8202 


100 


3 


10959.9356 


25794.5958 


18400 


9964.5902 


100 


4 


12570.2158 


24184.3126 


18400 


10367.1561 


100 


5 


13983.6886 


22770.8398 


18400 


10720.5172 


100 


6 


15200.3509 


21554.1774 


18400 


11024.6715 


100 


7 


16220.2029 


20534.3254 


18400 


11279.6153 


100 


8 


17043.2445 


19711.2838 


18400 


11485.3387 


100 


9 


17669.4757 


19085.0526 


18400 


11641.8108 


100 


10 


18098.8965 


18655.6319 


18400 


11748.8843 


100 


11 


18331.5069 


18423.0215 


18400 


11804.6769 



TABLE I: Energies of the bound states for the Woods-Saxon potentials for different values of n, I 
calculated using Eqs.(3.30), (3.33), (3.34) and (3.37). 
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